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Mott's metal-insulator transition at an interface due to band bending is studied by the density matrix 
renormalization group (DMRG) . We show that the result can be recovered by a simple modification 
of the conventional Poisson's equation approach used in semi-conductor heterojunctions. A novel 
mechanism of colossal electroresistance is proposed, which incorporates the hysteretic behavior of 
the transition in higher dimensions. 
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Strongly correlated electron system (SCES) - a group 
of materials in which the effect of Coulomb repulsion is 
large - is one of the major candidates on which the next- 
generation electronics may be built 0, [2, 0, 13 IE IE El E| • 
This strong hope put on to SCES electronics stems from 
the richness of its phase diagram 0; Various electronic 
and magnetic phases transitions are reported with high 
sensitivity to external conditions. Inside electric devices, 
the electrons behave collectively, so the high sensitivity 
of SCES may lead to drastic functionalities. 

Yet, SCES electronics faces a strong conceptual bar- 
rier to widespread acceptance and application, since it is 
essentially many-body where the useful concepts such as 
band structure are believed to fail. It is expected to be 
the case especially at the interfaces. In this Letter, we 
study interface Mott transition near an interface between 
a metal electrode and SCES in terms of the density ma- 
trix renormalization group (DMRG) method [Hj. sur- 
prisingly, we found that the conventional band bending 
picture based on Poisson's equation is valid with a small 
modification (eq. J3J). Namely, the conventional con- 
cepts in semiconductor devices are still useful and valid 
to design the SCES devices. As an application, we pro- 
pose a novel mechanism of the colossal electroresistance 
(CER), i.e., the large switching of resistance. 

Among many possibilities to industrial applications, 
perhaps CER in SCES hcterostructures is closest to real- 
ization, e.g., nonvolatile resistance random access mem- 
ories (RRAM) 0. The device consists of a film of 
perovskite manganite such as Pri-aCa^MnOs, <5 = 0.3 
(PCMO), which is a holed-doped Mott insulator, sand- 
wiched by two metallic electrodes. The current-voltage 
curve shows large hysteresis at room temperature, where 
the resistivity of the on and off states differ by a large 
factor. Although this effect has been explored in details 
IE HL IE j the understanding of the mechanism is still 
missing. Baikalov et al. pointed out from multileads re- 
sistance measurements, that the switching take place at 
the interface Q. Then, Sawa et al. reported that the 
CER behavior depends on the work-function of the elec- 
trode metal 0. They interpreted the I—V characteristics 



using a Schottky contact model (metal/p-type semicon- 
ductor) accompanied with an interface state. Another 
mechanism was proposed by Rozenberg et al. 

The mechanism of CER we propose here do not assume 
any interface states but attributes the large non-linearity 
of the I-V characteristics to interface Mott transition, 
where a layer of Mott insulator blocks the current. In 
dimensions higher than one, the transition is first order, 
and the width of the Mott insulator layer depends on 
how the voltage is changed. This explains the hysteretic 
behavior of the I-V characteristics. 

The Letter consists of two parts. We first study a 
one-dimensional model of a metal/SCES interface using 
DMRG combined with Poisson's equation. We show that 
interface Mott transition can be understood on the ba- 
sis of local equilibrium, i.e., electric state is determined 
by the local value of the potential. Then, we explore 
the hysteresis loop of the I-V characteristics on a phe- 
nomenological basis assuming the a hysteretic density- 
potential curve for a system in higher dimensions. Along 
with the mechanism of CER, another interesting conse- 
quence of interface Mott transition is proposed: A quan- 
tum well structure emerges spontaneously (Fig^a)) at 
the interface of a hole (electron) doped Mott insulator 
and an electrode with small (large) work-function. Our 
results may lead to fabrication of clean 2D metallic sys- 
tems analogous to the high electron mobility transistor 
in semiconductor physics. 

DMRG study of a ID interface — We start with the 
one-dimensional model of a metal/SCES interface on a 
lattice with total Hamiltonian 

-fftot = Hclc + #SCES + Hjnc, (1) 

where the electrode H c \ c = 

~t Ea E^l" 1 ( C ! + 1. C - + H - C ) and SCES 

#SCES = E<x Ef=L/2+l (4+la C ^ + H.cj + 

Eti/2+i {Un^Uii + ViUi) is connected by a junction 

H jnc = -*Ea i C \/ 2+ la C L/2a + H-c\ t is the hopping 
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element, L the total size of the system with the interface 
at the center, and U denotes the on-site Coulomb 
repulsion (for a related work, see 0). We take the 
units H = £o = 1. The potential V defined in the SCES 
region obeys ID Poisson's equation, whose discretised 
solution is 
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where m = (J2a c L Cicr ) ^ s ^ ne electron density and n + 
the positive background related to the hole doping ra- 
tio 6 by n+ =1 — 6, and a = (a : lattice const., e : 
dielectric constant), [a, 6] stands for the interval between 
a and b. -ffscESi without V, is identical to the Hamilto- 
nian of the one-dimensional Hubbard model, which ex- 
hibits a Mott transition at half-filling if U > [b|. If 
n =/= 1, the ground-state is metallic, a state known as the 
Tomonaga-Luttinger liquid (e.g. 

The DMRG calculation is performed as follows. In or- 
der to obtain the ground-state consistent with the poten- 
tial determined by Poissons's equation, we update the po- 
tential using eq.(0 at each step of the finite size method 
in the DMRG procedure ^3|- This is repeated more 
than 50 times to obtain total convergence. The typi- 
cal size of the block Hilbert space used here is m = 250. 
We fix the total number of electrons 'Y^=\ n i ~ Ln + , 
thus, the Fermi level of the electrode changes as electrons 
(holes) are injected to the SCES region when we vary Vqq. 
The interface is characterized by the work-function differ- 
ence Vd = Vl/2+i —Vl = (work-function of electrode) — 
(work- function of SCES). We note that Vx, achiving a 
given Vd depends on L. 

In Fig. ^ we plot the electron density and the potential 
in a mctal/SCES interface. Three typical solutions are 
plotted in a increasing order of Vd from (a) to (c). 

(a) Interface Mott transition: When the electrode's 

work-function is small enough, a quantum well is 
formed at the interface, i.e., two insulating (MI and 
BI) layers with width c?mi? <^Bii sandwiches a metal- 
lic region (M) with width c?m ■ The Vd dependences 
of these widths are plotted in Fig. [21 (b) . 

(b) Ohmic junction: The Fermi surface of the elec- 

trode and SCES balances and no barrier is formed. 

(c) Schottky barrier: A Schottky barrier is formed as 

in conventional metal/semiconductor interfaces. 

It is noted here that the qualitative features of the results 
are all captured well by the conventional band bending 
picture by replacing the valence (conduction) band by the 
lower (upper) Hubbard band as shown in Fig^ We also 
note that in the three cases, n{x) shows an oscillatory 
behavior in the electrode regime, which is the ID Friedel 
oscillation Sn(x) = cos(2fcpr + i]p)/r with kp the Fermi 
wave number, rjp a phase shift, and r the distance from 
the interface. 
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FIG. 1: Electron density n(x) (upper panel) and potential 
V(x) (lower panel) in a metal/SCES interface with (a) V x = 
-2.25, (b) Vx> = -1, and (c) V x = -0.5 (U/t = 4, n+ = 0.7, 
a = 0.043 and L — 100). In each panels, a description of each 
regions is given in the top where dark (white) regions are 
insulator (metal), e.g. band insulator (BI), metal (M), and 
Mott insulator (MI). The schematic picture in the left-lower 
space describes the bending of the upper (lower) Hubbard 
band UH (LH) near the interface. Ef is the Fermi energy. 



In FigEI(a), we plot the electron density n(x) against 
-(V(x) - Voo) in the SCES regime (x G [51,90]), where 
data near the boundary (x G [91, 100]) were omitted to 
avoid the boundary effect. The data falls onto a universal 
density-potential curve, which increases as the potential 
becomes deeper. In the middle of the curve, there is 
a plateau at half-filling whose width is A(U). The in- 
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FIG. 2: (a) Universal density-potential relation of the ID 
interface Mott transition. Electron density n(x), 51 < x < 90 
is plotted against —(V(x) — Voo). Voo is varied from Voo = 
-3.00 to Voo = with an interval of 0.25. W(U) is the band 
width and A({7) the width of the n = 1 plateau, (inset) 
The U dependence of A(U). This is determined from data 
with density n G [0.96, 1.04] (open circle) and n G [0.98, 1.02] 
(open box). The solid line is the Mott gap from Lieb-Wu's 
solution Hal , (b) The width of the Mott insulating (c(mi, 
n G [0.85,1.04]; blue), metallic (d M , n G [1.04,1.96]; red) 
and band insulating (dsi, n G [1.96, 2.0]; green) layers plotted 
against Vd. The symbols are DMRG results while the solid 
lines are the solutions of Poisson's equation (eqs. 0-©). 



terface Mott transition is an analog of the filling driven 
Mott transition There the filling n(fi) of a grand 

canonical system is studied while the chemical poten- 
tial /z is varied. The universal density-potential curve 
is an analogue of the n(fi) curve in the filling driven 



Mott transition, where —V(x), with a shift of the zero 
point, plays the role of /i. Indeed, when we compare 
the width of the plateau with the Lieb-Wu solution of 
the Mott gap of the one-dimensional Hubbard model 

A(t/) = W JT sin£^/u) d y El> the two coincides wel1 

as shown in the inset of Fig [3 (a). Thus, it can be said 
that the density-potential relation in the interface Mott 
transition follows the n(/i) curve in the bulk transition. 

Poisson's equation and local equilibrium approximation 
— Since the metal/SCES interface determines the trans- 
port properties of the device, the width of the layers c?mi, 
c^Mi and c?bi in Fig. ^ (a) is of practical interest. Here, 
we derive the width by solving the modified Poisson's 
equation 



d 2 V{x) 
dx 2 



-(n(V(x)) 

e 



(3) 



where we assume local equilibrium, that is, we assume 
that the electron density only depends on the local value 
of the potential. In the following, we evaluate eq. © 
to obtain the widths expressed solely by the potential 
difference Vb, hole doping ratio 5, band width W, and 
the Mott gap A. 

We adopt a simplified density-potential relation by lin- 
earizing the DMRG result Fig. [3(a); A constant com- 
pressibility — dn(V)/dV = k = 2/(W — A) is assumed for 
-(VXaO-Voo) € [{6-1)/k,S/k], [8/k+A,(1+5)/k+A], 
and n(V) = for -(V(x) - Vj) e [—oo, (5 - 1)/k], 
n(V) = 1 for -(V(x) - V^) G [S/k,S/k + A] and 
n(V) = 2 for -(V(x) - V^) G [(1 + S)/k+ A.oo]. We 
seek for a solution with a fixed density at the bulk SCES 
n(x) = n + = 1 — S, Mott insulator layer (MI) n(x) = 1 
and band insulator layer (BI) n(x) = 2, but varies in the 
metallic region (M). The width of MI layer is nonzero 
when Vd < —6/k 



1M1 



= ^2e(-V D - S/n)/eS 



(4) 



and saturates at Vd = —5/k — A. When the MI layer 
saturate, the metallic region appears whose width is 



^e/eKcosh -1 ((Vd + A)6/k + ^2A8/n v / (V D + A + 5/k) 2 - 2(S/k)(Vd + 6/k)X /{2A5/k - (6/k) 2 } , (5) 



which saturates when Vb < A. Finally, the width y/en/e (|) sinh {^/en/ed^Aj where d M here is obtained 

of the BI layer is nonzero when Vd < - A by substi t uting y D = _ I±£ _ A in cq. ©. In Fig. 0(b), 

rp , n-p x 2 oeTi t x\ TT\ T x\ i — it/ T~Al we P^ ^ ^ ne ^° dependence of the widths and compare 

d BI = 2 + V( 2) - -( +o){{ + )/k+Vd+ ) hem ith thc DMRG results and both agrees surpr i s _ 

1(1 + 6) 

ingly well. This agreement is highly nontrivial; In the 
insulating phase the localization length £ ~ W/A(U) is 



th F l = v / 2eA8/e and F 2 = F 1 cosh ( ^en/edyt J + f the order of few sites, but in the metallic phase it 
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should diverge. So, in the metallic phase the local ap- 
proximation is not a priori justified. However, our nu- 
merical calculation shows that it works remarkably well. 
Based on this success, we apply this local approximation 
to more generic cases below. 
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FIG. 3: (a) Density-potential relation near Mott's transi- 
tion in dimensions higher than one (schematic). A region 
(Vi < V < V2) exists where metallic and insulating phases 
coexist, (b) I-V characteristics of the metal/SCES interface 
when Vi < Va + Voo < V2 is satisfied, (c), (d) Pictures of the 
off state (c) and on state (d). 

Application to CER — In CER devices, the inter- 
face between the electrode and SCES is two dimensional. 
In systems with dimensions higher than one, the filling 
driven Mott transition is first order and hysteretic behav- 
ior takes place near the transition point [Tj, [l(| Il3 • We 
expect similar hysteresis in the density-potential relation 
of the interface Mott transition (Fig. El (a)). If this is 
the case, coexistence of the metallic state and the insu- 
lating state is possible. Then, the I-V characteristics of 
the device shows high non-linearity as well as a memory 
effect, i.e., RRAM behavior, when a layer of coexistence 



state exist in the SCES region. The condition for this 
to take place is Vb < — (6 — A6i)/k. For example, if 
Vi < Vd + Voo < V-2 is satisfied, the region nearest to the 
interface is a phase coexistence layer, as shown in Fig. 
13 (c) and (d). If this is the case, the junction may be 
insulating (c) or metallic (d). They correspond to the 
off and on states, respectively. We can switch between 
them by applying a voltage Va on the electrode forming 
a hysteresis loop in the I-V characteristics Fig. 01(b). 
We assume a tunneling form for the current 



/(Va) =T(V a )I (V a ), 



(7) 



where T(Va) is the tunneling probability of the bar- 
rier and Io(Va) the current at the metal/SCES junc- 
tion without any barrier. Iq(Va) reflects the details of 
the device and may be Ohmic, i.e. /o(Va) oc Va, or if 
space charge limited current is realized, -Zo(Va) oc (Va) 2 
|18l | . We assume that the tunneling probability decreases 
exponentially as the width of the MI layer grows, i.e., 
T(Va) = e ~ dMl{ y^^, where £ is the decay length (here 
the temperature dependence is neglected). Neglecting 
the jump of n in Fig. |3] (a), the width of the MI 
layer can be obtained by replacing — Vb — 5/ n in eq.(0} 
with Vi — Vd — Voo — Va , i = 1, 2 for the on and off 
state, respectively. Thus, a hysteresis loop is realized 
in the I-V characteristics as in Fig. [3] (b). If we de- 
fine the CER ratio by AR/R = (i? ff ~ Ron) /Ron, where 
Ron,oB is the resistivity of the on and off states, we ob- 
tain AR/R = eV^W-Vn-v^/eS/s _ 1. i n a Ti/PCMO 
based CER device, however, neither the on nor the off 
states show Ohmic I-V characteristics @. This can 
be also understood by our model with Vd + Voo < 
V\. In such cases, the CER ratio becomes AR/R = 

e ^2s(V 2 -Vn-V ao )/eS-^2e(V 1 -V D -V^)/eS)/^ _ ± ^ 

case, we can design a CER device with larger CER ratio 
by decreasing the doping ratio S and making the phase 
coexistence region wider. 

In summary, we have studied the interface Mott tran- 
sition by the DMRG method and by Poisson's equation 
combined with a local equilibrium ansatz. We proposed a 
novel mechanism of CER for materials with a first order 
metal-insulator transition. 
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